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Abstract 

A continuum limit of the non-Hermitian spin-1/2 chain, conjectured re- 
cently to belong to the universality class of the noisy Burgers or, equivalently, 
Kardar-Parisi-Zhang equation, is obtained and analyzed. The Galilean in- 
variance of the Burgers equation is explicitly realized in the operator algebra. 
In the quasi-classical limit we find nonlinear soliton excitations exhibiting the 
u; cx fc^ dispersion relation with dynamical exponent z = 3/2. 
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The Kardar-Parisi- Zhang (KPZ) equation plays an important role in modern non equi- 
librium statistical mechanics as a continuum description of growing interfaces [jl],|2[. In a 
co-moving frame the equation of motion for the relative height variable h{x,t) has the form 

dh/dt = uV^h + ^\{Vhf + r] , (1) 

where u is an effective diffusion coefficient, A a coupling constant character- 
izing the slope dependence of the growth velocity, and ri{x,t) a white noise, 
< ri{x,t)rj{x' ,t') >= A6{x — x')5{t — t'), representing fluctuations either in the drive or 
in the environment. 

The slope variable u = Vh satisfies the noisy Burgers equation 

du/dt - XuVu = z/V^M + Wt] . (2) 

The universality class represented by (|I]) and (0), besides driven interface and irrotational 
fluid dynamics, includes such important models as a driven lattice gas 0, as well as a 
directed polymer [Q, or a quantum particle 0], in random environment. These models serve 
as paradigms in the theories of driven and disordered systems. An additional attraction of 
the problem is its analytical tractability. A consistent, although uncontrolled, perturbative 
renormalization group flow has been constructed up to the two-loop level @J1,P] and appears 
to give qualitatively correct results. What emerges is a critical (massless) behavior, which at 
and below two spatial dimensions is always governed by a strong-coupling renormalization 
group fixed point. Two basic scaling dimensions describe the long- wave- length low- frequency 
behavior of the problem: i) the roughness exponent ( characterizing the height correlations 
in (|l]) and ii) the dynamical exponent z. The velocity field u or, equivalently, slope field V/t 
for the interface, has scaling dimension 1 — C- The pair correlation function scales as 

< u(xo, toMxo + X, to + t) >= a;"2(i-C) ^ (3) 

where the scaling function f{x) oc a;^^^^"^)/^ for large values of its argument. 

The case one spatial and one time dimension, = 1-1-1, is of particular interest. Although 
the perturbative renormalization group flows into the strong coupling regime, outside of its 
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range of validity, it was noticed in that both scahng dimensions follow explicitly from 
two exact scaling relations. The first one follows from the Galilean invariance of the original 
Burgers equation (generalized to A 7^ —1), i.e., u{x,t) — > u{x — Xuot,t) — uq, leaving all aver- 
ages unchanged. This property implies that A is a structural constant of the symmetry group 
and therefore remains invariant under a renormalization group transformation. Comparing 
the scaling properties of the two terms on the left hand side of we then obtain the first 
scaling relation, z = 2 — (. The other scaling relation is particular to the case d = 1 + 1 and 
follows from the fact that the simple Gaussian equal time stationary probability distribution, 
InP(u) oc — / u'^{x)dx, satisfies the Fokker-Planck equation for the probability distribution 
P{u,t), 



dP{u,t) 



dt 



dx 
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corresponding to (0) for any value of A. This property then implies ( = 1/2 (random walk 
behavior) and thus z = 3/2. 

The exact rational values of the scaling exponents found in d = 1 + 1 bear a superficial 
resemblance to the exponents encountered in two-dimensional classical critical problems. 



equivalent to ci = 1 + 1 relativistic quantum field theories |T^. Note, however, that the 
value 2; = 3/2 is manifestly non-relativistic, making extensions of the conformal field theory 
methods a highly nontrivial task. We might, however, still expect to gain new insight by 
studying the equivalent (l+l)-dimensional nonrelativistic quantum system by the methods 
developed for one-dimensional spin chains [|ll|]. Here spin waves and solitons in, e. g., the 
Heisenberg ferromagnet, are characterized hj uj = with z = 2 [Q, while activational 
quantum dynamics of the Ising chain in a random transverse field implies z = 00 [|T^ . 

Such a quantum spin chain approach has been proposed in and considered further 



in , on the basis of the equivalence between the master equation describing the evolution 
of a one- dimensional driven lattice gas or the equivalent lattice interface growth model and 
the non-Hermitian spin s = 1/2 Hamiltonian for spins. 



^ - 1 ^ ^ 1 
h = -J: Sr - - + Z6 ■ (5, X . (5) 

i=i ^ 

The mapping is achieved by identifying the eigenvalues of the z-projections of the spins, Sj, 
with the slope variable uj of the discrete interface model or the occupation numbers in the 
lattice gas representation. The interface dynamics is governed by two basic rates, r| and r^, 
corresponding to flipping "up" or "down" a kink (a pair of neighboring interface segments 
with opposite slopes); these flips map onto the spin exchange processes in The vector 
e is oriented along the 2;-axis and absorbing the mean rate + by rescaling time we 
have |e| = e = (r^ — rj)/(r| + r^), measuring the strength of the drive. For e = "up" 
and "down" flips are equally probable, and we obtain the Heisenberg ferromagnet, with 
the dynamical exponent z = 2, corresponding to the noisy linear diffusion equation or the 
Edwards- Wilkinson equation for interface dynamics [|l|,|l^. The asymmetric limit of only 
"up" or only "down" flips, e = ±1, has been solved exactly by the Bethe-ansatz method 



1^ , and the finite-size gap in the spectrum has been shown to scale with the length of the 
chain to the power 3/2, in agreement with z = 3/2 for the noisy Burgers universality class. 

In this Letter we obtain the continuum limit of the spin chain corresponding to the 
noisy Burgers universality class. We show that properly defined this limit exhibits Galilean 
symmetry as realized by the algebra of generators of the continuous global transformations. 



It is known that in the Heisenberg ferromagnet |12 , as well as in other XYZ spin-1/2 chains 
T^, quasi-classical quantization of the classical solutions of the continuum equations of 
motion correctly reproduces the low-energy sector of the Bethe-ansatz solution. Motivated 
by this fact we also study the quasi-classical limit. We find that the classical equations 
of motion have exact solitary wave solutions, which after quantization form a branch of 
elementary excitations with the cj oc k^^'^ dispersion relation. 

As a guide to the proper continuum limit we make use of the two exact relationships 
employed above in the derivation of the scaling dimensions ( and z. The implementation of 
the second one is readily carried out by noting that the stationary equal time probability 
distribution of the Burgers equation maps onto the ground state of the quantum problem. 



The ground state of the Heisenberg ferromagnet, |0), for e = is ahgned with maximum 
total spin sN and fully degenerate spin direction. In the following we quantize along the 
x-direction, i.e., (0|S'|'|0) = s, (015110) = 0, and (Ol^JlO) = 0, corresponding to a horizontal 
interface with zero slope. With this choice (cf. |T^) Sj is completely disordered and the 
correlation function {0\SjSj,\0) = s^djji. In the continuum limit ^j.yla <5(a;), where a 
is the lattice spacing, implying the rescaling S'^(x) ~ Sj/a}~^ with = 1/2. This result 
is also obtained by noting that a block spin S^{x) constructed from a random sequence of 
Si = ±1/2 yields a spin amplitude of the order of the square root of the block size. Finally, in 
the mapping of a master equation onto an equivalent quantum Hamiltonian the energy levels 
correspond to the relaxation rates. The ground state maps onto the stationary distribution, 
consequently, its energy must be equal to zero. It is easily seen that the aligned ground state 
|0) is also an eigenstate of H in (|^) with eigenvalue zero, owing to the cross product form 
of the last term. Consequently, the ground state structure and thus the scaling dimension 
of the spin density is not changed in the presence of the drive. 

The excitation spectrum, however, is altered since we expect the u oc k"^ dispersion law 
of the Heisenberg ferromagnet to change to oc k^^"^ in the presence of the drive term. 
Clearly, the Galilean invariance which permits the evaluation of the exponent z requires the 
continuum limit a — to be taken since discrete distances cannot be mixed with continuous 
time. 

With the ground state aligned in the x-direction the excitations above the ground state 
correspond to spin fluctuations in the y and 2;-directions. Noting that a block spin pointing 
in the x-direction has a size s ~ for a — * while the spins in the y and z directions scale 
as a~^/^, we infer that the continuum limit is equivalent to the large spin limit in the block 
spin representation. Using the quasi-classical harmonic oscillator representation for the spin 



operators, see e.g. [11], we have 
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S] = ^js^'\ = u,s^'^ . (6) 
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The oscillator position and momentum operators (pj and Uj, satisfy the commutator relation 
[(y5j,Mfc] = and correspond to the polar angle and the 2;-component of the spin, i.e., the 
action-angle representation of the Heisenberg spin [|1^]. In the continuum limit the field 
operators u{x) and (p{x) absorb the scale factor a^/^, u{x) = Uja'^/"^ and (y5(x) = (p^a'^^'^ , 
and obey the canonical commutation relation \(p[x\u[x')\ = i6{x — x'). Introducing the 
spin- wave stiffness J = a^^^, using periodic boundary conditions, and omitting constant 
terms, the Hamiltonian finally reads 



H = J dx 



(7) 



we note that the overall factor a^^"^ is consistent with the dynamical exponent z = 3/2. It 
is easily seen that the first term in (0) corresponds to the spin-wave approximation for the 
Heisenberg ferromagnet (the first term in (^); the second term representing the drive then 
corresponds to a spin-wave interaction. 

We can now show that the Hamiltonian (|^) is equivalent to the Fokker-Planck equa- 
tion (^) for the noisy Burgers equation (^. Introducing the conjugate operators u and 
i5/5u and performing the canonical transformation i5/5u = / A)^^"^ [ip{x) — iu{x)] and 
u = (A/z/)^/^'u(x) the operator C in dP{u,t) /dt = C(u)P{u,t), Eq. @ is equivalent to 



for J = u and e = |A(A/z/ 



1/2 



Noting that the time evolution operator for the master equation is exp(—tH) |T^, the 
equations of motion for and u follow from d0/dt = [H, 0] and du/dt = [Hju], 

dip ■j^'^'^ 2 "^'^ (8) 

dt dx"^ dx ' 

du ■j^'^'^ _|_ 2 " (9) 

dt dx"^ dx 

These field equations are invariant under the Galilean transformation 

u{x, t) — > u{x — 2euot, t) — uo , ^^^^ 
ip{x, t) — > ip{x — 2euQt, t) , 
in addition to being invariant under an arbitrary shift in 0. These two continuous symmetry 
relations correspond to the two generators of the full group of rotations of the Heisenberg 



spin in @ for e = 0. The oscillator representation corresponds to a local Abelian limit of 
the non- Abelian group 0(3); the spin wave Hamiltonian is invariant under arbitrary shifts 
in both canonical fields (p and u. In the presence of the spin wave interaction, the invariance 
under shifts in (p is preserved while the other symmetry is replaced by the Galilean space-time 
mixing ([TD|). 

These symmetries have a compact representation in terms of the algebra of the operators 
= J dxu and $ = / dxip generating rotations about the y and the 2;-axis, or, equiva- 
lently, shifts in ip and m, respectively. Adding the momentum operator P = J dxud/dip and 
H, generating space and time translations, respectively, we arrive at the operator algebra: 
[H, M'] = [H, P] = [P, $] = [P, IIP] = 0, and [l>, AP] = iL , where L is the length of the 
interface, while the Galilean invariance yields 

[H, $] = -2itP . (11) 

Thus e is indeed a structural constant of the algebra of symmetries. This justifies the 
renormalization leading to (|^) which were chosen so that e is invariant under changes in the 
microscopic scale a. The above commutation relation has a simple interpretation in the case 
of the elementary excitations, the lowest energy states carrying momentum k\ (|TT]) implies 
rotation of the global polar angle $ with constant angular velocity d^/dt = 2ek. 

We now attempt to find an elementary excitation with the to oc fc^^^ dispersion relation 
by a direct analysis of the nonlinear equations of motion (8) and (9). Such an analysis 
is most easily carried out in the classical limit followed by quasi-classical quantization (cf. 
||T^). We therefore first replace the operators and u in (|^) by the field variables (p and 
u satisfying the Poisson bracket {u{x), fiy)} = 6{x — y). The classical equations of motion 
du/dt = —i{H, u} and dip/dt = —i{H, ip} then have the same form as (8) and (9) with u and 
replaced by u and (p. For e = we find small amplitude spin waves about the ground state 
with quadratic dispersion leading to z = 2. However, for e 7^ we identify non-linear non- 
perturbative soliton solutions propagating with velocity v, i.e., replacing d/dt by —vd/dx 
we obtain —vdu/dx = —iJd'^ip/dx'^ + 2eudu/dx and —vdip/dx = iJd'^u/dx'^ + 2eudip/dx 
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which are easily solved by quadrature. 

However, for our purposes we note that by forming the quotient of the two equations 
and integrating once we obtain (du/dx)"^ + {dip/dxY = C. For solutions with vanishing 
derivatives at infinity C = 0, and we have 



By insertion we obtain one equation for the slope u which by quadrature yields —vu = 
±Jdu/dx + eu"^ + B Relating B to the boundary values u± = u{±oo) far away from the 
soliton we obtain du/dx = ±{e/J){u — u+){u — uJ) and the "mean amplitude- velocity" 
condition 



defining permanent profile kink or soliton solutions. The soliton shape is given by |(n — 
— = exp [±(e/J)(M+ — 'U_)(x — Xo)] with "center of mass" position xq. For 
du/dx = ±id(p/dx we have solitons (kinks) with positive and negative slopes, respectively. 
In order to interpret these solutions as excitations above the ground state, which classically 
is the uniform solution m = 0, we have to satisfy the boundary conditions m(±oo) = 0. We 
can achieve this by grouping kinks in pairs: for f > we group the = (0, —v/e) 

kink at smaller x with ('U_,-u_|_) = (— ti/e, 0) at larger x; for v < Q the = (0,t'/e) 

kink is followed by = (f/e, 0). The distance between the kinks in a pair is much 

larger than the intrinsic width J/v but much smaller than the length of the system. 

The energy momentum relationship for a pair is obtained by noting that the energy 
E is given by while the momentum, the generator of translation, is given by P = 
— / dxudip/ dx. Inserting ([T^ ) we obtain for a sing le kink, E = ±(1/3)03/2(^3 _ 

and 

P = ±(l/2)('u^ — v?l). Note that the two kinks in each pair moving with a constant velocity 
V correspond to the opposite choices of sign in (|12D ; as a result both energy and momentum 
are doubled and we obtain 



du/dx = ±id(f/dx . 



(12) 



u+ + u 



v/e , 



(13) 




(14) 
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for V > (<) 0. We note that the momentum is purely imaginary; this is a feature of the 
complex field equations (8) and (9). Eliminating the amplitude v /ewe obtain the dispersion 
law 

E oc at (±iP)^ . (15) 

In order to revert to a discussion of the quasi-classical limit of the spin chain the final step 
in our analysis amounts to a quasi-classical quantization of the above pair solutions. In 
the spirit of the Landau quasi-particle picture we envisage that we can label the low lying 
quantum states in terms of a dilute gas of pairs of non-linear kinks with dispersion law (|15|). 

A few comments are in order: a) The amplitude of a pair, f /e, diverges for e — implying 
that these are indeed non-perturbative solutions of the field equations characterizing the 
strong coupling renormalization group fixed point. This is also consistent with the fact that 
the coupling e does not appear in the dispersion law (p!5D, displaying the universality of the 
strong coupling regime, b) The kink solutions characterized by ( [T3| ) are manifestly invariant 
under the Galilei transformation (JT^). c) Inserting d(p/dx = idu/dx in the field equations 
(8) and (9) we obtain the deterministic noiseless Burgers equation. We conclude that the 
above kink solutions correspond to the well-known shock wave solutions ||^. In the interface 
representation a pair of kinks forms a step: a localized slope fluctuation propagating in the 
lateral direction. 

In conclusion, two main results for the theory of the noisy Burgers-KPZ universality 
class have been obtained in this Letter. First, although supported by an exact Bethe ansatz 
solution for a special value of the parameter e and by numerical results, the equivalence 
between the lattice model and the continuum Langevin equations (1) and (2) remained 



a conjecture |I^,|T5[. Here we have shown, using well-known methods in the quantum theory 
of spin chains |]1T|, how to implement the continuum limit of the lattice model (^). The 
rescaling leads to a continuum Hamiltonian (|^) which by a canonical transformation is 
reduced to the Liouville operator (^) of the Fokker-Planck equation corresponding to the 
noisy Burgers-KPZ equations. The field operators appearing in the continuum Hamiltonian 
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absorb the powers of the lattice constants making their dimensions equal to the scaling 
dimensions of the corresponding quantities of the noisy Burgers-KPZ universality class. The 
choice of the scaling dimensions is governed by the (generalized) Galilean invariance of the 
theory (|10]), which is explicitly realized in the commutation relation ( [Tl| ) of the symmetry 
algebra. 

Second, the space-time scaling t oc with the exponent z = 3/2 appears in the theory of 
the d = 1 + 1 noisy Burgers equation as the only one consistent with the Galilean symmetry of 
the problem [Q. However, little physical intuition exists concerning the physical realization 
of this dynamical scaling. Here we have made a step towards a better understanding of 
the physical nature of the problem by identifying explicit nonlinear soliton solutions to the 
classical limit of the nonrelativistic field theory. These solitons correspond to the shock 
wave solutions of the deterministic Burgers equation. Grouping shock waves of the opposite 
sign moving with the same velocity, we obtain a localized excitation with respect to a 
uniform solution, which in the KPZ representation corresponds to a step propagating along 
a growing interface. In the quasi-classical limit these soliton pairs form a branch of quantum 
excitations, exhibiting the nontrivial cu oc k^^"^ dispersion relation. 

Discussions with Jari Kinaret, Alan Luther, and Alexander Nersessyan are gratefully 
acknowledged. 
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